Detection of vacuum entanglement in a linear ion trap 
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We propose and study a method for detecting ground-state entanglement in a chain of trapped 
ions. We show that the entanglement between single ions or groups of ions can be swapped to the 
internal levels of two ions by sending laser pulses that couple the internal and motional degrees of 
freedom. This allows to entangle two ions without actually performing gate operations. A proof of 
principle of the effect can be realized with two trapped ions and is feasible with current technology. 



A remarkable phenomenon that appears naturally in 
quantum field theory, is that the ground state (vacuum) 
is entangled, and that observables in two separated re- 
gions can be entangled. Recent studies in quantum in- 
formation theory have taught us that entanglement is a 
physical property which can be exchanged between sys- 
tems or used in quantum processes such as quantum com- 
munication, teleportation and quantum cryptography 
This suggests that vacuum entanglement as well could 
be detected and used in quantum processes. There have 
been several studies of vacuum entanglement in field the- 
ory 0, 13 1 as weu as m other systems 0, IE @ , but none 
have proposed a way to observe vacuum entanglement in 
a realistic experiment. The main purpose of this Letter is 
to suggest a realistic physical implementation to observe 
this phenomenon. 

A gedanken-experiment that allows the observation 
of vacuum entanglement in field theory has been sug- 
gested Q, which utilizes two basic ingredients of rela- 
tivistic field theory and quantum information: the pres- 
ence of a causal structure, and the non-increase of entan- 
glement under local operations; operations performed at 
two causally disconnected regions do not increase the en- 
tanglement between these regions. Consider two atoms, 
A and B, which locally interact with the field and with 
one another through the long range field interaction. The 
interaction with the field can either entangle A and B via 
the exchange of propagating quanta, or by transporting 
vacuum entanglement into the atoms. Using the fields' 
built-in causal structure, one can eliminate the former 
unwanted process, by demanding that cT < L, where T 
is the interaction time and L the separation between the 
atoms. Vacuum entanglement can then be "swapped" to 
the atoms' internal levels, which can then be used for de- 
tecting vacuum entanglement. However this method re- 
quires precise control of the atom- field interaction, which, 
in the case of an electromagnetic field, renders the exper- 
iment highly unrealistic. Nevertheless, experimentalists 
in atomic physics are taming their systems at the quan- 
tum level and can test quantum effects with the required 
precision and control for this type of experiments. 

In the following, we propose and analyze the possi- 
bility of observing vacuum entanglement with trapped 
ions. We consider a system of trapped ions that are 
brought to equilibrium in a linear chain configuration. 
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FIG. 1: Scheme for detecting ground state entanglement: the 
entanglement between two groups of ions is swapped to the 
internal level of ions A and B by sending separate laser pulses 
that induce an interaction between the internal levels and the 
position of each ion. 



The ground state (vacuum) of the system is an entangled 
state of the different motional modes (phonons) of the 
chain, and manifests entanglement between single ions 
or distant groups of ions. In order to detect vacuum en- 
tanglement, we consider processes wherein the external 
motional degrees of freedom are mapped to the internal 
ions states, which are then used for entanglement detec- 
tion. The internal levels are well isolated, they can be 
temporarily coupled 'on demand' to the phonon modes 
by sending finite duration laser pulses, and finally can be 
measured with nearly perfect precision. In analogy with 
the field-theoretical case, the interaction must be limited 
to a duration shorter than the time it takes for pertur- 
bations to propagate between the two (probe) ions along 
the chain. We comment that in the case of ion chains, 
this process is interesting on its own, because one can 
entangle the internal levels of two ions without actually 
doing gates Q. The most spectacular manifestation of 
the idea would involve a chain with many ions. However 
a proof of principle can be attained with just two trapped 
ions, and is feasible using current technology. We shall 
study both cases and analyze the latter in detail. 

We now consider a system of N ions trapped in a lin- 
ear Paul trap at very low temperature @- The Hamil- 
tonian describing the ions' motion relative to their equi- 
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librium positions, and internal levels is Ho = \lo z {o z + 

Oj ) + ^2 VnO-ndn- Here lo z is the internal levels energy 
gap of the two relevant (probe) ions A and B, and v n are 
the phonon normal-mode frequencies, with correspond- 
ing creation (annihilation) operators a\ (a n ). Typically, 
lj z is in the optical region and v n ~ MHz. 

We shall begin by analyzing the simplest case with just 
two trapped ions. The vacuum state is then defined as 
the ground state of the normal modes of the system, i.e. 
a product state of the collective and breading modes |0 C ) 
and | Ob). In terms of the local single oscillator states 
\ri)A,B, the vacuum is an entangled two mode squeezed 
state 9] 

|vac> - |0 c )|0 b > = y/l-e-^J2e- 0n \n) A \n) B (1) 

n 

The local number states are the single ion energy eigen- 
states obtained when the displacement of the other ion is 
set to zero. We get e _/3 = a/(A - l/2)/(A+ 1/2), where 
A = (1 / '4)[^/uq I V\ + \Jv\j Vq] , and v$ and v\ are the fre- 
quencies of the collective and breathing modes. Since 
v$lv\ = y/3, we get A — 0.5189, and the von Neu- 
mann entanglement ^(j °f the squeezed state is E = 
(A + 1/2) log 2 (A + 1/2) - (A - 1/2) log 2 (A - 1/2) = 0.136 
e-bits. 

In order to transfer vacuum entanglement into the ion 
internal states, we use laser pulses to couple the internal 
and motional states of atoms A and B. Close to reso- 
nance, with LUiaser ~ Wzj the interaction terms for the 
k'th ion (in the Lamb-Dicke limit) is given by |Tlj 

HM=n{t){e-«*W+e»a<?>)x k , (2) 

where o~± are the raising and lowering operators, 4> is 
the laser phase, and Xk the displacement of the fc'th ion. 
Above we have applied the rotating wave approximation 
with respect to the internal levels but not (as is usually 
done) to position operators. That is because, the du- 
ration T of the laser pulses satisfies l/u> z -C T < 1/vq • 
The upper bound on T follows from the requirement that 
perturbations do not propagate between the ions during 
the interaction. 

Using the available interaction (2), we would like to 
swap the motional entanglement into the internal ion lev- 
els which are initially prepared in a non-entangled state 
I 1)a\ 1)b- The smallness of e"' 3 implies that the entan- 
glement arises mostly from the first two terms of • We 
therefore seek a procedure, of typical duration T -C l/voi 
that maps 

|vac)U)U)-| X )[||)U)+c^|T)|T)] (3) 

where x is the final state of the ions. The interaction then 
acts separately on each ion, and swaps the lowest two 
motional states |0) and |1) to the two ion internal states. 
In this 4x4 subspace the map is expected to approxi- 
mate the unitary swap transformation e™ '/ 4 ( cr = cr =+ cr y B) ^ 
where a x = |0)(1| + |1)(0| and a v = i\l){0\ - £|0)(1| act 
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FIG. 2: Histogram of the final density matrix of the internal 
levels of ions A and B. The entanglement of formation of this 
state accounts for 97% of the computed initial ground state 
motional entanglement. 

on the number states. (We have ignored a <j z o~ z term 
which in our case adds a trivial phase since the initial 
state is a | J,)| [) state). We note that a x and 5 y can 
be approximated by x and p respectively. We therefore 
expect that the swap will require coupling of the inter- 
nal levels with both position and momentum, since the 
truncation of the x operator yields the a x operator and 
truncation of momentum yields the a y operator. Based 
on this intuition we proceed with the following construc- 
tion. We consider the following sequence of unitary oper- 
ations, U s = V(a 1 )W(p 1 )V(a 2 )W((3 2 )---V{a n )W((3 n ), 
where 

V(a) = e ia(T * x W{p) = e l(3 ^ p , (4) 

to be performed on each ion separately by sending a se- 
quence of laser pulses. The V(a) evolution, can be ob- 
tained by sending a laser pulse of duration T and phase 
4> = 0, such that T <C 1/vq an d / £l(t)dt = a. In or- 
der to generate a W{(3) evolution, we set the laser phase 
to = 7r/2, and allow the system evolve freely for short 
time interval dt — t in between a pair of pulses. Denoting 
V'(/9) = exp(i(3a y x), we obtain: 

VU(-P')V^ (P') = e -*(/3'M* + ^)+0(, 2 T 2 ) e */3W 

= g-i^pr+lr/J') + {u\ 2 ) (5) 

where we have used the approximation x(t) — x(0) + 
p(0)r/m + 0{v 2 t 2 ). (Alternatively, in the Schrodingcr 
picture we notice that e~ lHphonont V shifts the kinetic 
term as p 2 -> (p + (3'o- y ) 2 = p 2 + 2/3'a y p + (3 12 .) Taking 
the limit v 2 t 2 <C 1, and maintaining (3 = P'r/m = 0(1) 
we obtain the required effective coupling to p. There- 
fore the sequence of n pairs of V(a)W((3) pulses, can be 
generated by 3rt ordinary pulses with n free evolution 
intermediate intervals of total duration dt = nr. By op- 
timizing the entanglement of formation [T^. Ep(cti, Pi), 
over the free parameters on and Pi the transformation 
©, where x is the final motional state of the ions, 
can be generated with high efficiency. This transfor- 
mation swaps the first two terms in Eq. to the 
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FIG. 3: Entanglement (in e-bits) between complementary 
symmetric groups of ions as a function of the total ion num- 
ber. 

ion's internal level states. After a sequence of three VW 
pulses, the entanglement of formation of the final inter- 
nal level state contains 97% of the computed ground- 
state entanglement. The optimal sequence is in this case 
V{0M)W(0.38)V(0.50)W(0.39)V(0.53)W{0.16). (With 
two pulses we get at most 93%). Expressed in the rele- 
vant 4x4 subspace, this unitary operation has indeed a 
a structure which closely resembles a swap. Testing the 
purity of the final state pab, we find trp\ B = 0.997. The 
final density matrix of the internal levels is depicted in 
Fig. |21 The measurement precision for the density matrix 
in recent experiments is about 1% 0, and hence suffi- 
cient for observing the entanglement of the final state. 

We comment that an alternative to the above approach 
could be to separate the two ions by moving them apart. 
This effectively "turns off" the interaction between them, 
allowing for longer duration of the detection process. It is 
then easier to generate the desired swap using continuous 
on-resonance laser pulses 01- I n order not to affect the 
entanglement between the ions, the separation has to be 
fast compared with the propagation time scale 1/Vo- 
This can be done by increasing the potential between 
the ions. The possibility of changing the local potential 
and moving ions, without effecting the internal ion states, 
has been recently demonstrated experimentally [TH Il6| . 

We next turn to the more general case of entanglement 
in a chain with N ions, and consider first the entangle- 
ment properties of the system. Since the vacuum is a 
pure Gaussian state, the Schmidt decomposition for two 
complementary parts of the chain can be expressed as a 
direct product of squeezed states [ij . The von- Neumann 
bi-partite entanglement can then be obtained by sum- 
ming the entanglement contribution of each squeezed 
state (Fig. 0). 

The state of two sub-groups, A and B, each consist- 
ing of tia and tib ions, separated by l s ions, is de- 
scribed by a reduced Gaussian density matrix p^ B = 
t r i^A s(l vac )( vac D- The entanglement between the 
groups (Fig. 0} can be characterized by the Negativ- 
ity fl7l. Ilq. It vanishes for separation larger then one. 
However as the group size increases, it persists for larger 
separations. 

We shall examine the possibility of detecting the en- 
tanglement by coupling to the internal levels of two ions, 
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FIG. 4: Logarithmic Negativity between two groups consist- 
ing of 1,3, and 5 ions, as a function of their separation in a 
chain of 20 ions. 

A and B, for a time duration T. We can check that as 
long as T < 1/Vo, the interaction with respect to two 
complementary parts of the chain can be regarded as lo- 
cal. This is seen by noticing that the evolution operator, 
[/(T), can be factorized in the interaction picture as 

U = U A ® U B ® e- 1 ' 2 I dtdt'f(t,t')a AB ^ (g) 

where U k act on A or B, o A b = n^e^o^ + e^cr <>) ), 
and f(t— t') = [xA{t),XB (t')]. The last term above, 
involving ctab, is a unitary operator that can increase 
entanglement "non-locally" . However, as can be seen 
in Fig. [SI the non-commutativity described by f(t — 
t') vanishes rapidly, and for sufficiently short interaction 
time, or large enough spatial separation, this non-causal 
effect is suppressed. 

We begin with the initial ground-state |vac)| |)| |), 
and proceed to evaluate the reduced state pab(T) = 
tr[U(T)\I)(I\W(T)} perturbatively. Assuming that the 
intensity of the laser pulses is sufficiently weak, we ex- 
pand U(T) in a power series, and to lowest order in f2 we 
obtain 

(\\x AB \\ 2 -{o\x AB ) \ 

||Ba|| 2 (e b \e a ) 

(E A \E B ) \\E B f 

~(x AB \o) 1-||Ba|| 2 -||Bb|| 2 / 

(7) 

Here \E A ) = Xi|vac), \X AB ) = X A X B \vaa) and X k 
| dttt(t)e iSt x k (t), (k = A, B), and 5 is the detuning. Us- 
ing the Peres-Horodecki separability criterion, it then fol- 
lows that pab{T) is entangled iff N(pab) ~ \ {0\Xab) \ ~ 
\\Ea\\\\Eb\\ > 0, where N(pab) is the negativity. This 
condition can be understood physically as the require- 
ment that the virtual off-shell single phonon exchange 
process (described by Xab), is sufficiently large to over- 
come the decoherence effects due to local phonon emis- 
sion (described by the Ea } b terms). 

To verify that the above condition amounts to a de- 
tection of vacuum entanglement, rather than a direct in- 
teraction due to the non-local correction in Eq. (3), we 
have repeated the computation, using the same initial 
state |7), but with a modified truncated evolution. In 
the latter truncated case, we "disconnected" the chain 
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FIG. 5: a) Classical propagation of a perturbation originat- 
ing at the center of the chain, b) The commutation relation 
between the displacement operators of the n'th ion and the 
central ion in a chain of 80 trapped ions, at different time 
slices. 




FIG. 6: The ratio r\ for ions in a chain of TV = 20 ions as a 
function of the detuning 8. a) The ion probes are located at 
I = 6, 15 and T = 0.8 (in unites wherein i>o = 1). b) I = 10, 11 
and T — 0.05. The range r\ > 1 signifies entanglement. 

by eliminating the interaction between ions at the differ- 
ent halves, n > N/2 and n < N/2, of the chain. This can 
be easily achieved by replacing the potential term in the 
free phonon Hamiltonian by ^XiGijXk — > ^XiGj^Xk 
where G T = Ga © Gb is block-diagonal. This trun- 



cated evolution, does not change the entanglement be- 
tween the two halves of the chain since an exact separa- 
bility holds in Eq. (3), (i.e. f(t - t') = 0). The ratio 
r\ = \(p\X AB )\/\\E A \\\\E B \\ is plotted Fig. as a func- 
tion of the detuning S, for N = 20 ions. In Fig. [Ui the 
probe-ions are situated at sites I a — 6 and Ib = 15, and 
we can see that A and B become entangled (77 > 1) in a 
range of frequencies. Since small violations of causality 
are expected in the non-truncated model, one could have 
anticipated that the non-truncated case should give rise 
to more entanglement. On the contrary, we see that it is 
the truncated case which gives rise to more entanglement. 
To understand this consider first Fig. |B|d, in which near- 
est neighbor ions Ia,b — 10, 11 have been used as probes. 
We find that the truncated and non-truncated models 
precisely agree for sufficiently small T <C 1/Vo- In this 
case, since there is pre-existing local entanglement be- 
tween the close ions, the probe can detect entanglement 
in an arbitrarily short time, and the truncated interaction 
has here no effect because evolution is not required. On 
the other hand, in the case of separated ions, propaga- 
tion effects "communicate" between the probes and ions 
closer to the center of the trap which carry the most en- 
tanglement. This suggests that the larger entanglement 
in the truncated model is due to perfect wave reflection 
at the boundary between the regions. 

In conclusion, we have proposed an efficient method 
for detecting vacuum entanglement by mapping motional 
states of trapped ions or groups of ions to the ions' inter- 
nal levels. It is remarkable that this phenomenon, which 
can be considered as purely fundamental, may also help 
to realize quantum information tasks. Further investiga- 
tions will determine whether vacuum entanglement can 
be used to entangle internal degrees of freedom in a fast 
way, and to produce spin squeezing, which would be of 
practical interest in, for example, atomic clocks ^J. 

AR and BR acknowledge the support by ISF, grant 
62/01-1. JIC was supported by EU projects, the DFG, 
and the Bayerischen Staatsregierung. 



[1] M. Nielsen, I. Chuang, Quantum Computation and 

Quantum Information., Cambridge Univ. Press, 2000. 
[2] S. J. Summers and R. F. Werne r, Phys. Lett. A 110 , 257 

(1985); R. Verch, R. F. Werner, quant-ph/0403089 
[3] B. Reznik, Fo und. Phys. 33, 167 (2003); 

quant-ph/0008006 B. Reznik, A. Retzker, J. Silman, 

quant-ph/0310058 
[4] T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 

032110, (2002); G. Vidal et. al. Phys. Rev. Lett 90, 

227902 (2002). 

[5] K. Audenaraert et. al, Phys. Rev. A 66, 042327 (2002); 

A. Botero, B. Reznik, Phys. Rev. A 70, 052329 (2004). 
[6] U.V. Poulsen, T. Meyer, M. Lewenstein, 

cond-mat/0405550 
[7] I. J. Cirac, P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 
[8] D. Leibfried et. al, Rev. Mod. Phys. 75, 281 (2003). 



[9] A. Botero, B. Reznik, Phys. Rev. A 67, 052311 (2003). 

[10] A. Peres, Quantum Theory: Concepts and Methods 
(Kluwer, Dordrecht, 1995). 

[11] J. I. Cirac, L. M. Duan, P. Zoller, quant-ph/0405030 

[12] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[13] M. A. Rowe et al, Nature 409 ,791 (2001). 

[14] B. Kraus, J. I. Cirac, Phys. Rev. Lett. 92, 013602 (2004). 

[15] M. A. Rowe et al, Quant. Inf. Comp. 4 257 (2002). 

[16] M. D. Barrett et al, Nature 429 , 737 (2004). 

[17] G. Vidal, R. F. Werner, Phys. Rev. A 65, 032314 (2002). 

[18] The logarithmic negativity measures the degree of which 
the partial transposed density matrix (which is a trans- 
pose of the density matrix with respect to one of the 
subsystems) fails to be positive, moreover it bounds the 
distillable entanglement. 

[19] D. J. Wineland, et al, Phys. Rev. A 50, 67 (1994). 



